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Abstract 

The formulas for the temperature dependences of the non-strange and strange 
quark condensates are derived by taking into account the contribution of the 
massive resonances. Critical temperature of the chiral symmetry restoration 
transition is established to be 190 MeV if only meson resonances are consid- 
ered, and 175 MeV if both meson and baryon resonances are taken into account. 
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For massless quarks, the QCD Hamiltonian is invariant under chiral transformations. 
The phenomenology of the strong interactions is consistent with this property pro- 
vided that chiral symmetry is spontaneously broken. Lattice gauge calculations 
indicate that, at low temperatures, the scalar currents qq = {uu, dd, ss} develop 
nonzero expectation values (qq)r which are referred to as quark condensates, and 
represent order parameter of spontaneously broken chiral symmetry. Chiral pertur- 
bation theory 0, |3|, £| shows that (qq)r melts as the temperature increases. One 
generally expects that if the temperature reaches a critical value, T c , chiral symmetry 
is restored. Lattice calculations || confirm this expectation but show that the strange 
quark condensate (ss) does not melt in the chiral restoration transition, supporting 
the description of this phase transition in terms of the behavior of the SU(2) x SU (2) 
linear a-model initiated by Wilczek et al. 0. Presumably, deconfinement transition 
which liberates color, takes place at temperature of the order of T c . At low enough 
temperatures, chiral symmetry provides constraints on the temperature dependences 
of several physical quantities. For example, for the pion gas without interactions, one 
finds 1 

(qq) T _ 1 Nj-l T* f_J^_ = ^ _ JV)7* f^) =1+ 1 T 2 



(qq) 3N f 4/2' f n 6 4/*' m n 6JV/4/2' 1 

where f w is the pion decay constant, ~ 93 MeV, n n the screening pion mass, and Nf 
the number of massless flavors. In a real world, Nf = 2, so that the first relation of 
(1) may be rewritten as 

(99>t , T 2 x _( T 



(qq)o 8/2 V260 MeV 

If one goes beyond the low temperature limit, predictions start to become model 
dependent. 

In this paper we calculate the correction to the formula for the non-strange quark 
condensate, Eq. (2), given by the contribution of the higher mass hadrons, and ob- 
tain a similar relation for the strange quark condensate. It is well known that the 
correct thermodynamic description of hot hadronic matter requires consideration of 
higher mass excited states, the resonances, whose contribution becomes essential at 
temperatures ~ 0(100 MeV) [0, [8|, |9|. The method for taking into account these res- 
onances was suggested by Belenky and Landau [0 as considering unstable particles 
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on an equal footing with the stable ones in the thermodynamic quantities; e.g., the 
formulas for the pressure and energy density in a resonance gas readQ 

p = E* = £*^E<±ir'*^. (3) 

i i Z/l r=l ' 

P = E/ ) " Pi=T^-p t , (4) 

i 

where +1 (—1) corresponds to the Bose-Einstein (Fermi-Dirac) statistics, and gi are 
the corresponding degeneracies (J and / are spin and isospin, respectively), 

(2Jj + l)(2/j + 1) for non — strange mesons 
4(2Jj + 1) for strange (K) mesons 

2(2J i + l)(2J i + l) forbaryons, 

and since chiral symmetry suppresses the interactions of low energy Goldstone bosons 
both among themselves and with massive hadrons, the gas can be approximately 
described collection of free particles. 

To derive an explicit relation for the value of the quark condensate as a function 
of temperature, we first note that the operator qq occurs in the quark mass term of 
the Hamiltonian, 

H = H + J d 3 x m i 99> ( 5 ) 

q=u,d,s 

where H is the Hamiltonian of massless (chirally symmetric) QCD. The thermal 
expectation value of qq represents, therefore, the response of the partition function 
to a change in the quark mass, 

1 d\nZ 

{(iq)T = -v~d^> (6) 

where V is the four- dimensional euclidean volume (= V^'/T). At sufficiently large 
volume, In Z — > V(p — po), where p is the pressure and po the vacuum energy density. 
In the large volume limit, Eq. (6) therefore takes on the form 

(a) = dpo 9p = (- Q ) dp (7) 
T dm q dm q dm q 

Since p depends on m q only through the masses of the particles, one obtains from 
(3), (7), through d/dx (x 2 K 2 (x)) = -x 2 i^i(x), 

^^o+E^^Biir 1 ^^ 1 - (8) 



1 Since the temperatures we are dealing with are much less than the nucleon mass, we may treat 
all the particles as bosons, and distinguish fermions only by the factor 7/8 in the expression for the 
particle degeneracy. We also neglect chemical potential for simplicity. 
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1) Non-strange quark condensate. 
One uses the lowest order relation ITT 



m lfl = - m (q(i)o, (9) 
where (qq)o = (uu + dd) and m = 1/2 {m u + m<f), and obtains from (8), 
/- \ /- \ (i T 2 9i rrhmdrrii ^ +1 mi /rm,\\ 

1) Strange quark condensate. 
Now one uses the lowest order relation for the K meson (fx being the kaon decay 
constant, ~ 114 MeV), 

m 2 K f K = -1/2 (m + m s )(ss + 1/2 qq) = -m s (ss) , (11) 

since m s » m, and (ss) ~ (uu) ~ (dd)o [|T^j, and obtains from (8), 

Gerber and Leutwyler M have established the T-dependence of the non-strange quark 
condensate by direct calculation of the sum in Eq. (10) for known resonances, using 
the estimate drrii/dm ~ N(, the number of valence quarks of type u and d, and m ~ 7 
MeV. The results of Gerber and Leutwyler predict chiral symmetry to be restored 
at T c ~ 200 MeV if only massive states are taken into account, and with the pion 
contribution in addition to that of the massive states, T c ~ 190 MeV for nonzero 
masses of the u- and <i-quarks, and 170 MeV for zero masses of the latter. In this 
paper we wish to obtain analytic expressions for both the T-dependent non-strange 
and strange quark condensates, and compare the results with those of Gerber and 
Leutwyler. 

In the following we shall restrict ourselves to the meson resonances alone. Baryon 
resonances may be treated in a way similar to that described in this paper, with the 
introduction of two chemical potentials, for both conserved net baryon number and 
strangeness. We shall dwell briefly on this point at the end of the paper. 

To calculate the derivatives dmi/dm q in Eqs. (10), (12), we first note that the 
expansion of m% and m 2 K in the powers of m u , ma, m s starts with terms which are 
linear in the quark masses, 

ml = 2mB, m\ = (m + m s )B . (13) 

For higher mass meson resonances, we expect that similar relations hold (with C 
being constant within a given meson nonet): 

m 2 Q n ~ m\ = 2mB + C, m\ /2 = (m + m s )B + C, ml,~2m s B + C, (14) 
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where mi, mi/2, rn i, m » are the masses of the isovector, isospinor and two isoscalar 
states, respectively, and belongs to a mostly octet. Indeed, for the vector meson 
nonet these relations may be easily obtained from the constancy of the differences of 



the squared masses of the corresponding spin-triplet and spin-singlet states |13 

AM 2 

resulting inR 



M\ A S 1 



M 2 ( 1 S I 



m 2 p -ml 



0.57 GeV' 



m 2 K , - m\ 



0.55 GeV 2 



and two Gell-Mann-Okubo mass formulas, the standard one 14 



m\ + 3m\ = 4m 2 



and an extra relation 15, 16 



2 i 2 



1/2) 



m 8 + m = 2 m l/2 



/2' 



(15) 
(16) 

(17) 



with mo and m§ being the masses of the isoscalar octet and singlet states, respectively, 



which for an almost ideally mixed nonet reduces to [15, [16 



m 2 // ~ m 2 



l) 



m 2 / ~ 



2m\ 



/2 



m 1 . 



;is) 



Therefore, it follows from (13), (15), (18) that 



2 2 

m, , ~ m„ 



2mB + C, m 2 K * = (m + m s )B + C, m\ ~ 2m s B + C, 



with C ~ 0.56 GeV 2 . It was shown by Balazs JT^ that m 2 — m 2 = l/2a , with a 
being a universal Regge slope, a' ~ 0.85 GeV -2 , in agreement with (15), so that 



m 2 u) ~m 2 = 2mB + l/2a, m l K * = (m+m s )B + l/2a , ~ 2m s B + l/2a . (19) 

For higher spin nonets, since the corresponding states with equal isospin and alter- 
nating parity lie on the linear Regge trajectories, one has, e.g., 



m 



U2 



III 



m , 



rn 



m 2 p + l/a = 2mB + 3/2a, 
m 2 K *+l/a= (m + m s )B + 3/2a , 
m 2 + l/ct' ~ 2mB + 3/2a\ 
ml + l/a ~ 2m s B + 3/2a , 



(20) 



2 Also, m|,» - m|j = 0.55 GeV 2 , m 2 D , - m^^ = 0.58 GeV 2 , 
exception of the cc states for which "i 2 /^ — — 0.70 GeV 2 . 



= 0.55 GeV 2 , with the 
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m 2 p3 = m 2 a2 + l/a = 2mB + 5/2a, 

m K * = m K * + l/a = (m + m s )B + 5/2a , 

m w 3 = m / 2 + V a — 2mB + 5/2a, 

m 2 = m 2 , + l/«' ~ 2m s B + 5/2a, etc., (21) 

J 2 



and also 



m\ x = + 1/a = 2m5 + 1/a , 

m ifi = n^lc + l/ 02 = ( m + m s)B + 1/a', etc. (22) 

Thus, we consider the relations (14) as granted by both the Gell-Mann-Okubo mass 
formula for a close-to-ideally mixed nonet and the Regge phenomenology. It then 
follows from (13), (14) that 



dm » dmi mB m\ 



m " ~ m 



m 



dm dm mi 2m 

dmi/2 mB m 2 

dm 2m\/2 ^iii-i 
dm Q > 
dm 



m—^- ~ 0, (23) 



and 



dm n " dmi 
m 9 „ ~ m 



m. s 



<9m s dm s 
dmi/2 m s B m s m\ ^ m~ K 



dm s 2mx/2 2m 1/2 (m + m s ) 2m 1/2 



dm^_ _ r?vB _ m s m\ _ m\_ 
dm s m Q > m <(m + m s ) m ' 

since m s » m. Thus, we find that the expressions to be inserted into Eqs. (10) and 
(12) are, respectively, 

dmt ml dmi m\ . . 

m— — = , m s — — = , (25) 

dm 2mi dm s mi 

and out of 9 isospin degrees of freedom of a nonet, 6 contribute to the formula for 
the non-strange condensate: 3 isovector, 2 isospinor and 1 isoscalar which belongs to 
a mostly singlet, and 3 to the formula for the strange condensate: 2 isospinor and 1 
isoscalar which belongs a mostly octet .0 We shall make another simplification in Eqs. 
(10), (12), viz., approximate the Bose-Einstein statistics by the Maxwell-Boltzmann 



3 In fact, each of the uu, dd, ss condensates gains the contribution of 3 isospin degrees of freedom 
of a nonet. 
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one, taking account of the sum over r through the factor 7r 2 /6 = ((2) = X^l/r 2 , 
which is the asymptotic form of this sum for T » rrii. Then, the formulas (10) and 
(12) will finally reduce to 



(W>T = (qq)o 1 " T^E^^i h£ , (26) 




= <«>«> U - y^E^ 1 (t)J ' (27) 

where the primes indicate the summation over 6 and 3 isospin degrees of freedom 
of a nonet, respectively. Now it is seen in Eq. (26) that if one restricts himself to 
the massless pions alone, g n = 3, one obtains the formula (2). This formula sets 
the temperature scale for the non-strange condensate ~ 250 MeV. As seen in Eqs. 
(26), (27), the massive states accelerate the melting of the condensates, so that the 
temperature scale for, at least, the non-strange condensate is expected to be much 
narrower. 

To calculate the sums in Eqs. (26), (27), we resort to the notion of a resonance 
spectrum which is introduced in order to substitute the summation over individual 
particle species by the integration over the mass in the expressions for thermodynamic 
quantities, so that, e.g., Eqs. (3), (4) may be rewritten (in the Maxwell-Boltzmann 
approximation for the particle statistics) as 

rm,2 m?T 2 / m\ 

p = J dm r(m)p(m), p(m) = ~l^r K 2 [r^J » ( 28 ) 

p= dmr(m)p(m), p(m)=T — p( m ), (29) 

J mi dl 

and the resonance spectrum r(m) is normalized as 

J mi 



dm r(m) = ( 30 ) 



i 



where mi and m2 are the masses of the lightest and heaviest species, respectively, 
entering the formulas (3), (4). 



In both the statistical bootstrap model (18], [L9j and the dual resonance model (20 
a resonance spectrum takes on the form 

r(m) ~ m a e m/T °, (31) 

where a and To are constants. The treatment of a hadronic resonance gas by means 
of the spectrum (31) leads to a singularity in the thermodynamic functions at T = Tq 



18, [OJ and, in particular, to an infinite number of the effective degrees of freedom in 
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the hadron phase, thus hindering a transition to the quark-gluon phase. Moreover, 



as shown by Fowler and Weiner [21|, an exponential mass spectrum of the form (31) 
is incompatible with the existence of the quark-gluon phase: in order that a phase 
transition from the hadron phase to the quark-gluon phase be possible, the hadronic 
spectrum cannot grow with m faster than a power. 

In our previous work 2^] we considered a model for a transition from a phase 
of strongly interacting hadron constituents, described by a manifestly covariant rel- 
ativistic statistical mechanics which turned out to be a reliable framework in the 
description of realistic physical systems [p3fl , to the hadron phase described by a 



resonance spectrum, Eqs. (28), (29). An example of such a transition may be a rel- 
ativistic high temperature Bose-Einstein condensation studied by the authors in ref. 
23, which corresponds, in the way suggested by Haber and Weldon E3, to sponta- 



neous flavor symmetry breakdown, SU(3)p — » SU(2)j x [7(l)y, upon which hadronic 
multiplets are formed, with the masses obeying the Gell-Mann-Okubo formulas WA 



m 



a + bY + c 



Y 2 



(32) 



here / and Y are the isospin and hypercharge, respectively, I is 2 for mesons and 1 
bor baryons, and a, b, c are independent of I and Y but, in general, depend on (p, q), 
where (p, q) is any irreducible representation of SU(3). Then only the assumption on 
the overall degeneracy being conserved during the transition is required to lead to 
the unique form of a resonance spectrum in the hadron phase: 



r(m) = Cm, C = const. 



(33) 



Zhirov and Shuryak [^| have found the same result on phenomenological grounds. 
As shown in ref. PSfl , the spectrum (33), used in the formulas (28), (29) (with the 
upper limit of integration infinity), leads to the equation of state p, p ~ T 6 , p = p/5, 
called by Shuryak the "realistic" equation of state for hot hadronic matter 0, which 
has some experimental support. Zhirov and Shuryak p6| have calculated the velocity 



of sound, c s = dp j dp = c 2 (T), with p and p defined in Eqs. (3), (4), and found that 
c 2 (T) at first increases with T very quickly and then saturates at the value of c 2 s ~ 1/3 
if only the pions are taken into account, and at c 2 ~ 1/5 if resonances up to M ~ 1.7 
GeV are included. 

We have checked the coincidence of the results given by the linear spectrum (33) 
with those obtained directly from Eq. (3) for the actual hadronic species with the 
corresponding degeneracies, for all well-established hadronic multiplets, both mesonic 
and baryonic, and found it excellent [p^j . Therefore, the theoreticalconclusion that 
a linear spectrum is the actual spectrum in the description of individual hadronic 
multiplets finds its experimental confirmation as well. In our recent papers [15] , [Uj 
we have shown that a linear spectrum of an individual nonet is consistent with the 
Gell-Mann-Okubo mass formula (16) (in fact, this formula may be derived with the 
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help of a linear spectrum |T5|]), and leads to an extra relation for the masses of the 



isoscalar states, Eq. (17), which was checked in ref. [16| and shown to hold with an 
accuracy of up to ~3% for all well-established nonets. In ref. [IE we have general- 
ized a linear spectrum to the case of four quark flavors and derived the corresponding 
Gell-Mann-Okubo mass formula for an SU(4) meson hexadecuplet, in good agree- 



ment with the experimentally established masses of the charmed mesons. In ref. |27] 
we have applied a linear spectrum to the problem of establishing the correct qq as- 
signment for the problematic meson nonets, like the scalar, axial-vector and tensor 
ones, and separating out non-gg mesons. In this paper we shall apply a resonance 
spectrum to the derivation of the formulas for the temperature dependences of the 
quark condensates. 



It was shown in ref. |22] that the actual resonance spectrum does not depend on 



the dimensionality of spacetime.0 Therefore, the sums in Eqs. (26), (27) which are, in 
fact, related to the expression for the pressure of a free gas in 1+1 dimensions [29]], as 
calculated for individual nonets, may be substituted by the integration over m with 



a linear spectrum. The normalization constant C was established in ref. ||22[| : for a 
nonet, one has 9 isospin degrees of freedom lying in the interval (m » ~ mi, m '). 
Therefore, Eq. (30) gives 



and hence 



f m o' 
C / dm m 



C = — Z~2~ = — ~ 27 GeV" 2 , (34) 



mjy — m 2 ,, A 



where the difference A = m 2 , — m 2 „ is determined by a distance between the parallel 
Regge trajectories for the oj and resonances, which are described by the straight 
lines J = 0.59 + 0.84M 2 and J = 0.04 + 0.84M 2 , respectively, so that A ~ 0.65 GeV 2 . 
We note further that the meson nonets may be arranged in the pairs of nonets which 
have equal parity but different spins (which differ by 2), e.g.,0 

1 3 Pq jpc = 0++; ao (i320), /o(1300), / (1525), K*(U3Q), 

1 3p 2 jpc = 2++; a2 (i320), / 2 (1270), ^(1525), K* 2 (U30), 

1 3 D 1 J PC = 1— , p(1700), w(1600), A"*(1680), (no candidate), 
1 3^ jpc = 3 -- p3 (i700), w 3 (1600), 03(1850), K*(1780), 

3 1 S J PC = 0"+, vr(1770), 7/(1760), X(1830), (no 77' candidate), 
1 l D 2 J PC =2-+, vr 2 (1670), K 2 (1770), (no 77,77' candidates), 
and occupy the mass interval of an individual nonet but have 18 isospin degrees of 

This is another argument against the Hagedorn spectrum, since the exponent a in Eq. (31) 
depends explicitly on the dimensionality of spacetime (it is related to the number of transverse 
dimensions of a string theory p8j ). 

5 For the scalar meson nonet, we use the qq assignment suggested by the authors in ref. Jl6|, |27|. 
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freedom in this interval, i.e., twice as much as that for an individual nonet. Moreover, 
as the temperature gets closer to the critical one of chiral symmetry restoration, 
we expect the chiral partners (the states with equal isospin but different parity) 
have equal masses and form parity doublets. The work of DeTar and Kogut [jBOfl 
shows convincingly that the "screening masses" of chiral partners are different below 
and become equal above a common T c . This work was carried out for four quark 
flavors. Similar results were obtained for two flavors by Gottlieb et al. |3l| . In these 
calculations, the chiral partners were (vr,cr), (p, a±) and (iV(|+), iV(|— )). Thus, we 
expect the correct density of states per unit mass interval to be twice as much as that 
for an individual nonet, and hence, the correct normalization constant is 

C ~ 54 GeV" 2 . (35) 

In the case we are considering here, one has 6 isospin degrees of freedom in the interval 
m i/2 — m o" = ^V^> * n y i ew °f (17), contributing to the formula for the non-strange 
condensate, and 3 in the interval m 2 / — m\i 2 = A/2 contributing to the formula for 
the strange condensate; i.e., a lower half of a linear spectrum of a nonet contributes 
to (qq)r, while an upper half to (ss)t- The corresponding normalization constants 
are 

19 24 4 

C = — — = — = -C ~ 70 GeV" 2 , (36) 
qq A/2 A3 V ; 

C Ss = -V = ^ = -C ~ 35 GeV" 2 . (37) 
A/2 A3 v ; 

Once the mass spectrum of a nonet (with a given fixed spin) is established to be 
linear, one may take into account different nonets with different spins in Eqs. (3), (4). 
As shown in ref. [p2| , since the particle spin is related to its mass, J« ~ a mf, a 
being a universal Regge slope, the spin degeneracy turns out to be proportional to 
the mass squared, and the account for different nonets results in the following mass 
spectrum, 

T{m) = C'm\ d = 2a C ~ 90 GeV" 4 , (38) 

which is the actual resonance spectrum of hadronic matter and leads to the equation 
of state ^4 

p,p~T 8 , p = p/7. (39) 

Bebie et al. have calculated the ratio p/p directly from Eqs. (3), (4), with all known 
hadron resonances with the masses up to 2 GeV taken into account, and found that 
the curve p/p first decreases very quickly and then saturates at the value of p/p ~ 7, 
as read off from Fig. 1 of ref. 0, in agreement with (39). 

In order to show that the obtained normalization constant is correct, we note that 
the number of states with the masses up to M, given by the mass spectrum (38), is 

C' M 4 

N(M) = —— ~ 22.5 (M, GeV) 4 . (40) 
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For, e.g., M = 1.25 Eq. (40) gives 



7V(1.25) ~ 55. (41) 

The masses up to 1.25 GeV have the members of the pseudoscalar and vector meson 
nonets, and the /ii (1170), &i(1235) and ai(1260) mesons, the mass of the latter was 
indicated by the recent Particle Data Group as 1.23 GeV p2|. We do not include 



the scalar mesons ao(980) and /o(980) which seem to be non-gg objects |27j], but may 
include the /o(1300) meson which has the mass lying in the interval 1 — 1.5 GeV, 
according to the recent Particle Data Group. Thus, we have 9+27+1+9+9+1=56 
actual mesonic species having the masses up to 1.25 GeV, in excellent agreement with 
(41). 

For M — 1.7 GeV, Eq. (40) gives 

iV(1.7) ~ 188. (42) 

As seen in the Meson Summary Table [Q, the masses up to 1.7 GeV have the 
members of the following nonets: 1 %, 1 3 Si, 1 1 P U 1 3 P , 1 3 ^i, 1 3 ^2, 2 1 S , 2 3 Si. 
Therefore, one has 

(20 spin states) x (9 isospin states) = 180 states, (43) 

in good agreement with the result (42) given by a cubic spectrum. 
For M = 2 GeV, Eq. (40) gives 

iV(2) ~ 360. (44) 



The masses up to 2 GeV have the members of all the nonets indicated in except 
for the 1 3 i<4 and 2 3 P2 nonets. In this case, one has 

(41 spin states) x (9 isospin states) = 369 states, (45) 

again in good agreement with the result (44) given by a cubic spectrum. Thus, 
we consider the cubic spectrum (38) as granted by the actual experimental meson 
spectrum. 

Similarly to the case of an individual nonet, one finds the normalization constants 

C' qq = 2aC qq ~ 120 GeV" 4 , C'- ss = 2aC Ss ~ 60 GeV" 4 , (46) 

and the relations (10), (12) finally take on the forms, respectively, 

,_ . A T 2 120 GeV" 4 T 2 r , 3 rn fm\\ 
t = <«>o J dm m *- Kl (-) j , (47) 

<«>t - <«>o (l - 60 %^ JdmA, (^)) , (48) 
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where we have separated out the contribution of the pions which may well be treated 
as massless at temperatures ~ 150 MeV, and taken into account the remaining 
hadronic species by the resonance spectrum (38) with the normalization constants 
given in (46). We have found that the contribution of the kaons (as well as r\ and 
j] ), if separated out as well in Eqs. (47), (48), is comletely negligible at these tem- 
peratures. We include these states together with the other mesonic states in the 
integrals of Eqs. (47), (48), which, therefore, have the lower limit of integration ~ 0.5 
GeV. With respect to the remaining integration in (47), (48), we note that the main 
contribution to integrals of this type is given by the mass region in which m » T; 
therefore, one may extend the upper limit of integration to infinity and neglect the 



lower limit, and obtain, through the formula [35 



2"-V- 1 r 



dx x^K v [ax) ■ 

(qq)t - (qq)o 

(ss) T ^ (ss) 



1 + // + zA /! + // 



T 



T 



260 MeV/ V220 MeV 

rp \ 6 

1 



235 MeV 



(49) 
(50) 



Temperature dependences of the condensates, as given by (49), (50), are shown in 
Fig. 1. The critical temperature of the chiral restoration transition is T c ~ 190 MeV, 
in agrement with the result of Gerber and Leutwyler obtained for the pions and 
massive states with the nonzero masses of u- and ci-quarks. For the strange quark 
condensate, the temperature at which it is melted out is ~ 235 MeV, as seen directly 
in (50). For T ~ 190 MeV, Eq. (50) gives (ss) 190 ~ 0.72 (ss) , i.e., only 28% of 
the strange quark condensate melts by the chiral restoration transition. It is seen 
in Eq. (49) that the share of the pions in the overall reduction of the non-strange 
condensate from its initial value at T = to zero at T c is about 53%, the remaining 
47% is the contribution of the massive states. We note also that Eq. (49) may be 
well approximated by 



(qq) T ~ (qq) 



1 



T 



190 MeV 



(51) 



The formulas (49), (50) have been obtained for the meson resonances alone. There 
is no difficulty of principle to consider the baryon resonances in a similar way, with 
the inclusion of two chemical potentials, for both conserved net baryon number and 
strangeness. As we have checked in ref. ||22|| , a mass spectrum of the SU(3) baryon 



multiplets is linear, as well as for the meson nonets, although to establish its corre- 
spondence to the Gell-Mann-Okubo formulas is more difficult than for a meson nonet, 
since these formulas are linear in mass for baryons (more detailed discussion is given 
in ref. ||22||). Recent result of Kutasov and Seiberg |3H| shows that the numbers of 
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bosonic and fermionic states in a non-supersymmetric tachyon-free string theory must 
approach each other as increasingly massive states are included. The experimental 
hadronic mass spectrum shows that in the mass range ~ 1.2 — 1.7 GeV, the number 
of baryon states nearly keeps pace with that of meson states [0 (and, therefore, is 
well described by the same cubic spectrum as for the mesons, Eq. (38)). Above ~ 1.7 
GeV, the number of the observed baryons begins to outstrip that of the mesons, and 
then greatly surpasses the latter at higher energies (indicating, therefore, that the 
baryon resonance spectrum grows faster than (38) in this mass region, since the cubic 
spectrum (38) describes the meson resonances well, up to, at least, 2 GeV, as we 
have seen in Eqs. (41)- (45)). The explanation of this behavior of the experimental 
resonance spectrum was found by Cudell and Dienes in a naive hadron-scale string 
picture |3"Sfl: the ratio of the numbers of the baryon and meson states should, in fact, 



oscillate around unity, with the mesons favored first, then baryons, then mesons again, 
etc. Keeping in mind this picture, we may assume that the "in-average" baryon reso- 
nance spectrum has the same form, Eq. (38), as the meson resonance one. It is then 
possible to estimate the contribution of the baryon resonances to the temperature 
dependences of the quark condensates by assuming that the formulas (26), (27) hold 
for the baryon resonances, as well as for meson ones, and that out of 27 degrees of 
freedom of SU(3) baryon octet, nonet and decuplet, 2/3, i.e., 18, contribute to (qq), 
and 1/3, i.e., 9, to (ss). If one now neglects, for simplicity, the baryon number and 
strangeness chemical potentials (i.e., considers the case of both zero net baryon num- 
ber and strangeness), and takes into account the baryon resonances along with the 
meson ones by the mass spectrum (38), one will obtain the same formulas, (47), (48), 
but with the factors in front of integrals which are twice as much as those in the case 
of the meson resonances alone. These formulas will further reduce to the relations 



(qq) T ~ (qq) { 



1 



T 



260 MeV 



T 



195 MeV 



(52) 



(ss) T ~ (ss) 



T 



210 MeV 



(53) 



shown in Fig. 2. One sees that now the critical temperature of the chiral symmetry 
restoration transition is T c ~ 175 MeV, while the strange quark condensate melts out 
at ~ 210 MeV, and about 34% of the latter melts by the chiral restoration transition. 
Now the pions' contribution to the overall reduction of the value of the non-strange 
condensate is about 45%, the remaining 55% is the share of the massive states. 



Concluding remarks 

We have derived the formulas for the temperature dependences of the quark conden- 
sates by taking into account the contribution of the massive states parametrized by a 
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resonance spectrum which has been established by the authors in previous papers. In 
the case of the meson resonances alone, our results agree those obtained previously 
by Gerber and Leutwyler, and suggest the chiral symmetry restoration transition to 
occur at the critical temperature T c ~ 190 MeV. With the baryon resonances taking 
into account along with the meson ones, the critical temperature turns out to be 
T c ~ 175 MeV, which is closer to the currently adopted value 140 — 150 MeV estab- 
lished by lattice gauge calculations |39|. In either case, only ~ 30% of the strange 



quark condensate melts by the chiral restoration transition, in agreement with lattice 
calculations [||]. The share of the pions in the overall reduction of the non-strange 
condensate from its initial value at T = to zero at T c is ~ 50% in either case, the 
remaining ~ 50% is the contribution of the massive states. 



Acknowledgements 

One of us (L.B.) wish to thank E.V. Shuryak for very valuable discussions on hadronic 
resonance spectrum. 



References 

[1] For a review, see S. Gottlieb, Nucl. Phys. B (Proc. Suppl.) 20 (1991) 247, 
B. Petersson, Nucl. Phys. A 525 (1991) 237c, 
F. Karsch, Nucl. Phys. A 590 (1995) 367c 

[2] J. Gasser and H. Leutwyler, Phys. Lett. B 184 (1987) 83 

[3] J. Gasser and H. Leutwyler, Phys. Lett. B 188 (1987) 477 

[4] P. Gerber and H. Leutwyler, Nucl. Phys. B 321 (1989) 387 

[5] J.B. Kogut, D. Sinclair and K.C. Wang, Phys. Lett. B 263 (1991) 101 

[6] R.D. Pisarski and F. Wilczek, Phys. Rev. D 29 (1984) 338 
F. Wilczek, Int. J. Mod. Phys. A 7 (1992) 3911 
K. Rajagopal and F. Wilczek, Preprint PUPT-1347 (1992) 

[7] E.V. Shuryak, The QCD Vacuum, Hadrons and the Superdense Matter, (World 
Scientific, Singapore, 1988) 

[8] C.R. Alcock, G.M. Fuller and G.J. Mathews, Astrophys. J. 320 (1987) 439 

[9] H. Bebie, P. Gerber, J.L. Goity and H. Leutwyler, Nucl. Phys. B 378 (1992) 95 



14 



[10] S.Z. Belenky and L.D. Landau, Sov. Phys. Uspekhi 56 (1955) 309; Nuovo Cim. 
Suppl. 3 (1956) 15 

[11] M. Gell-Mann, R.J. Oakes and B. Renner, Phys. Rev. 175 (1968) 2195 

[12] M. Shifman, A. Vainshtein and V. Zakharov, Nucl. Phys. B 147 (1979) 385,488 
L.J. Reinders, H. Rubinstein and S. Yazaki, Phys. Reo. 127 (1985) 1 

[13] W. Lucha, F.F. Schdberl and D. Gromes, Phys. Rep. 200 (1991) 127 

[14] S. Okubo, Prog. Theor. Phys. 27 (1962) 949, 28 (1962) 24 

M. Gell-Mann and Y. Ne'eman, The Eightfold Way, (Benjamin, N.Y., 1964) 

[15] L. Burakovsky and L.P. Horwitz, Gell-Mann-Okubo Mass Formula for SU(4) 
Hexadecuplet, to be published 

[16] L. Burakovsky and L.P. Horwitz, Mass Spectrum of a Meson Nonet is Linear, to 
be published 

[17] L.A.P. Balazs, Phys. Lett. B 120 (1983) 426 

[18] R. Hagedorn, Nuovo Cim. 35 (1965) 216; Nuovo Cim. Suppl. 6 (1968) 169, 311 

[19] S. Frautschi, Phys. Rev. D 3 (1971) 2821 

[20] S. Fubini and G. Veneziano, Nuovo Cim. A 64 (1969) 811 

[21] G.N. Fowler and R.M. Weiner, Phys. Lett. B 89 (1980) 394 

[22] L. Burakovsky, L.P. Horwitz and W.C. Schieve, Hadronic Resonance Spectrum: 
Power vs. Exponential Law. Experimental Evidence; to be published 

[23] L.P. Horwitz and C. Piron, Helv. Phys. Acta 46 (1973) 316 

L.P. Horwitz, W.C. Schieve and C. Piron, Ann. Phys. (N.Y.) 137 (1981) 306 
L.P. Horwitz, S. Shashoua and W.C. Schieve, Physica A 161 (1989) 300 
L. Burakovsky, Manifestly Covariant Relativistic Statistical Mechanics as a 
Framework for Description of Realistic Physical Systems, Ph.D. thesis (Tel-Aviv 
University, 1995), unpublished; L. Burakovsky, L.P. Horwitz and W.C. Schieve, 
Mass - Proper Time Uncertainty Relation in a Manifestly Covariant Relativistic 
Statistical Mechanics, to be published 

[24] L. Burakovsky, L.P. Horwitz and W.C. Schieve, A New Relativistic High Tem- 
perature Bose-Einstein Condensation; Phys. Rev. D, in press 

[25] H.A. Haber and H.E. Weldon, Phys. Rev. D 25 (1982) 502 (1994) 4725 

[26] O.V. Zhirov and E.V. Shuryak, Sov. J. Nucl. Phys. 21 (1975) 443 



15 



[27] L. Burakovsky and L.P. Horwitz, Hadronic Resonance Spectrum May Help in 
Resolution of Meson Nonet Enigmas, to be published 

[28] K.R. Dienes and J.-R. Cudell, Phys. Rev. Lett. 72 (1994) 187 

[29] H.A. Haber and H.E. Weldon, J. Math. Phys. 23 (1982) 1852 

L. Burakovsky, On Relativistic Statistical Mechanics, Thermodynamics and Cos- 
mology in 1+D Dimensions, Tel-Aviv Univ. Preprint TAUP-2187-94, unpublished 

[30] C. DeTar and J. Kogut, Phys. Rev. Lett. 59 (1987) 399; Phys. Rev. D 36 (1987) 
2828 

[31] S. Gottlieb et al, Phys. Rev. Lett. 59 (1987) 2247 
[32] Particle Data Group, Phys. Rev. D 50 (1994) 1429 
[33] Particle Data Group, Phys. Rev. D 50 (1994) 1320 

[34] L. Burakovsky and L.P. Horwitz, On the Thermodynamics of Hot Hadronic 
Matter, in preparation 

[35] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, (Aca- 
demic Press, New York, 1980), p. 684, formula 16 

[36] D. Kutasov and N. Seiberg, Nucl. Phys. B 358 (1991) 600 

[37] P.G.O. Freund and J.L. Rosner, Phys. Rev. Lett. 69 (1992) 765 

[38] J.-R. Cudell and K.R. Dienes, Phys. Rev. Lett. 69 (1992) 1324 

[39] See, e.g., C. Bernard et al, Phys. Rev. D 45 (1992) 3854 



16 



FIGURE CAPTIONS 



Fig. 1. The ratio (qqjr/ '{qq)o as a function of temperature for a) non-strange, b) 
strange quark condensate, in the case of the meson resonances alone, as given in Eqs. 
(49),(50). 



Fig. 2. The same as Fig. 1 but in the case of both the meson and baryon reso- 
nances, as given in Eqs. (52), (53). 
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